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^'. Abstract 

(D . 

'"^ ' We consider models with N U{1) gauge fields A" N Kalb- Ramond fields 

-B" , an arbitrary bare action and a fixed UV cutoff A. Under mild assump- 

^ . tions these can be obtained as effective low energy theories of SU(N-|-1) Yang 

H I Mills theories in the maximal abelian gauge. For a large class of bare actions 

they can be solved in the large N limit and exhibit confinement. The confining 
phase is characterized by an approximate "low energy" vector gauge symme- 
try under which the Kalb-Ramond fields -B" transform. The same symmetry 
allows for a duality transformation showing that magnetic monopoles have 
condensed. The models allow for various mechanisms of confinement, de- 
pending on which sources for A^ or 5" are switched on, but the area law for 
the Wilson loop is obtained in any case. 



LPT Orsay 01-69 
June 2001 



lUnite Mixte de Recherche - CNRS - UMR 8627 



1 Introduction 

Monopole condensation is widely believed to be the mechanism responsible for 
confinement in Yang-Mills (YM) theories [|l], |^ (see [^ for a review). The picture of 
a confining vacuum as a dual superconductor gives rise to models for the low energy 
behaviour of YM theories which are based on (the dual of) the abelian Higgs model 
1^. Albeit quite successful phenomenologically P, 0, it is practically impossible to 
derive these models in a systematic way from YM theories: By construction (i.e. 
the assumption of duality) the gauge and scalar fields of the dual models are related 
non-locally to the YM gauge fields, at least off-shell, hence local dual models can 
only capture the semi-classical features of low energy YM theories. 

Generalizing a formalism proposed by Julia and Toulouse Quevedo and Tru- 
genberger [Q have presented quadratic actions for antisymmetric tensor fields, which 
describe the condensation of various topological defects in various space-time di- 
mensions d. If apphed to monopole condensation in d = 4, this approach suggests a 
(local) action for a Kalb-Ramond field B^^, [Q. 

There are two ways to understand the usefulness of 2-form fields B^^, for the 
description of monopole condensation in d = 4: First, the dual abelian Higgs model 
in the broken phase necessarily has to contain a pseudoscalar Goldstone boson. 
Reversing the duality transformation in d = 4, the dual of a pseudoscalar is a 2- 
form field B^^. Second, in the presence of a monopole the abelian Bianchi identity 
9[fj,Fi,p] = is violated on the world-line x{t), and the field strength F^jy cannot 
be written in the form F^u = d[f^A,^] in x(t). If monopoles have condensed in the 
vacuum, the abelian Bianchi identity is violated everywhere, and the field F^^, is 
nowhere of the form d^^A,^]. In this situation it is useful to introduce an auxiliary 
field Bf^,y for F^j, |]lOl in some analogy to the field strength formulation of YM theories 
|TT[| (it suffices, however, to introduce auxiliary fields S^,^ for the components of F^^ 



associated to the Nc — I U{1) subgroups of SU{Nc), see below). 

All this does not imply that a 2-form field -B^j^ appears in the physical spectrum 
of a theory with condensed monopoles: The equations of motion for -B^^, as derived 
from the full effective action, can well be algebraic, i.e. the propagator of B^^, has 
not necessarily a simple pole for finite negative p^ (in the Minkowskian regime). The 
role of the B^^ dependent terms in the full effective action is then just to parametrize 
in a compact way a subclass of higher derivative interactions among the fundamental 
gauge fields A^, but the modes in B^,^ which are not of the form d[fj_A^j correspond to 



topologically non-trivial gauge field configurations A^ associated to a non-vanishing 
monopole density. 

Monopole condensation is a purely abelian phenomenon even in the context of 
YM theories. Abelian dominance becomes particularly clear |T^ in the maximal 
abelian gauge (MAG) ||13[: the N^ — 1 gauge fields of SU{Nc) are divided into the 
N = Nc — I U{1) gauge fields A^, n = 1 . . . N, associated to the generators of 
the ?7(1)^ Heisenberg subalgebra, and A^^ — N^ "charged" gauge fields W^ (which 
carry at least one of the A^ U{1) charges) associated to the off-diagonal generators. 
The continuum MAG is characterized by the gauge condition (Df^Wf^)"' = 0, where 
Dfj^ are the f/(l)^-covariant derivatives. Lattice results indicate that the W^^ gauge 
fields are massive in the MAG, with M\y ~ 1 GeV (for SU{2)) |]I^. A mechanism 



for Wfj^ mass generation based on ghost condensates has been proposed in [1^ . 

Combining the above considerations one is led to conclude that the natural low 
energy degrees of freedom of confining gauge theories are N U{1) gauge fields A"^, 
and A^ Kalb-Ramond fields -B" which serve to describe monopole condensation. 
(Here "degree of freedom" does not necessarily signify an asymptotic physical state, 
but possibly just an auxiliary field with algebraic equations of motion.) 

The purpose of the present paper is the solution of a field theory with the cor- 
responding field content, with a bare action as general as possible, in the large A^ 
limit. The idea is to treat this theory as an effective low energy theory valid below 
some scale A ~ Myy (with, e.g., A ~ 1 GeV). Consequently we will supplement the 
model with an UV cutoff A, and allow for arbitrary non-renormalizable interactions 
~ A^^ (where p is determined by power counting) in the bare action. The precise 
numerical value of A will not concern us here. 

As we will find below, these models indeed describe confinement provided the 
parameters of the bare action satisfy some inequality. In the confining phase a 
duality transformation to an abelian Higgs model can be specified (although the 
full dual action is still non-local), which shows that one describes indeed monopole 
condensation. Remarkably the effective action exhibits an approximate vector-like 
"gauge symmetry" (a slight exaggeration, since it is broken by derivative terms) 
which allows to "gauge away" the low momentum modes of the U{1) gauge fields A" 
They are swallowed by the fields -B" in analogy with the Higgs-Kibble phenomenon 
with an additional Lorentz index. This symmetry of the low energy effective action 
coincides with the one in for the quadratic terms of the action. Here, however, 
it holds for part of the interactions (to arbitrary powers in the fields) as well and. 



simultaneously, renders a more general duality transformation possible. 

The conventional criterium for confinement is the area law of the expectation 
value of the Wilson loop. In an abelian theory the Stokes theorem allows to express 
the Wilson loop in terms of a source localized on a 'Wilson surface" bounded by 
the Wilson loop. In our model such a source can be coupled either to F^i, (denoted 
by J-^), to Bfj_i^ (denoted by J^) or both of them. As we will discuss in detail, in the 
confining phase of the model the area law is obtained in any case. We will show, on 
the other hand, that only a particular linear combination of the sources J^ and J^ 
is consistent, once the model is coupled to external quantum fields: Fluctuations of 
such quantum fields would generate an infinite contribution to the action unless the 
long-range parts of the correlators are cancelled. The remaining short-range parts of 
the correlators coincide with various models for the expectation value of Yang-Mills 
field strength correlators P, |l6i . 

A full-fledged computation of the string tension would require the solution of 
a coupled set of equations of motion (in analogy to approaches based on the dual 
abelian Higgs model [^) which will not be attempted here. In any case the result 
will be "non-universal" in the sense that it depends on details of the bare action 
and the UV cutoff. 

In this paper we focus rather on the infrared regime of a conflning gauge theory 
and the remarkable fact that, given the above fleld content representing the feature 
of "abelian dominance" and a non-perturbative computational scheme like the 1/A^ 
expansion, such theories can actually be solved in spite of confinement. 

The subsequent sections of the paper are organized as follows: In section 2 we will 
discuss more precisely, how the model can be obtained from SU{Nc) YM theories, 
specify its field content and parametrize its bare action. Section 3 is devoted to its 
solution at large N. Since the visibility of the confining phase necessitates a careful 
treatment of the limit of an infrared cutoff k"^ tending to zero, we will first review 
the bosonic 0{N) model at large A^ in the broken phase: Already in this simple 
model a convex effective potential, and hence a consistent treatment of the broken 
phase, requires the introduction of /c^ 7^ and a discussion of non-analyticities for 
k"^ -^ 0. Subsequently the emergence of a confining phase of the A^ — B^i, - model is 
seen to share similar non-analyticities with the 0{N) model, which originate again 
from the convexity of the effective action (and not from the large A^ limit). 

In section 4 we will discuss various properties of the confining phase: the possi- 
bility of a duality transformation, the emergence of a vector-like approximate gauge 



symmetry, the effective propagators of the A^ — B^i, - system, the "response" of 
the system to external sources like the Wilson loop and the unsuccessful search for 
bound states (hence "glueballs" will only emerge once the charged gluons of the 
SU{Nc) YM theory is taken into account). In section 5 we compare the large N 
limit of our models to the large Nc limit of SU{Nc) YM (which turns out to be not 
the same), discuss generalizations of the bare action considered in section 3 as well 
as the notion of universality, and conclude with an outlook. A short summary of 



the subsequent results has already appeared in ^5 



2 The A^ - B^^ - model 

Before giving a detailed description of our model, we will motivate it starting 
with the pure SU{Nc) YM theory in the MAG. As noted in the introduction, one 
splits the gauge fields into the N = Nc — I U{1) gauge fields A^ and the charged 
gauge fields PV^. The Euclidean partition functions reads 

1 f 
-^ VA VWVcVc e-^y«-s,/-s,h ^2.1) 

where the YM action can be decomposed as 

Sym = J d'x {f;^f;, + f;,f;4 (2.2) 

and M' is a normalization. 

The gauge fixing terms are of the form 

with the [/(l)^-covariant derivative 

D^W; = d,W; + gf:,A-^Wl . (2.4) 

The form of the ghost interactions in Sgh does not play any role subsequently. 

Now we introduce Kalb-Ramond fields 5^^, as collective fields for the N field 
strengths F^^^ [including their non-abelian terms ~ g f^^ ^hi^u]) ^^ the partition 
function ( p.l|) . This amounts to the multiplication of the integrand of ( p.l|) by 

l = ^jVBe-'- , 5B = /rf^a:{l(F;,-i?;,)'} . (2.5) 



The new version of the partition function reads 

— fvAVB VWVcVc e-^^^'-^sf-^a'^-^^ (2.6) 

with TV = M'Mb- 

Next, as stated in the introduction, we assume that the VW path integral is 
"infrared save" : Lattice resuhs show that the ly^ gauge fields in the MAG aquire a 
finite mass M\y, in agreement with the phenomenon of infrared abelian dominance in 
this gauge [Q. A mechanism for W!^ mass generation which involves the four ghost 



interactions in Sgh in the MAG has been proposed in |T^. The same mechanism (bi- 
ghost condensation) renders the "charged" ghosts c", c" massive. Now we assume 



that the VW and VcDc path integrals in (|2.6| ) are performed. However, in order not 



to generate uncontrollable UV divergences, it is advisable to integrate simultaneously 
over the "high momentum" modes of A" and 5" : These modes are needed in order 
to cancel UV divergences which would violate the Slavnov- Taylor identities. 
Let us thus split the measures VA and VB as 

VA = VA^M^-VA^M^ , VB = VB^M^-VByM^ . (2.7) 

The indices <Mw and >Mw denote modes with momenta in the corresponding 
ranges. Now we rewrite the partition function (|2.6|) as 



A/'< 
with 



y J VA^M^ VB^M^ e-S,ar.{A,B,Mw)-S,f (2. 



K 



s,^,.^(A,B,Mw)-s,f = _sJvA^M^ VB^M^ VW VcVc 



p—SYM—Sgf—Sgf^—Sl, 

(2.9) 

Sgf in (|2.8| ) and on the left hand side of (|2^ ) corresponds to the first (abelian) term 
in ( p. 31) . The partition function ( p.8|) corresponds to the class of models which are 
the subject of the present paper. It represents the "effective low energy theory" of 
SU{Nc) YM under the only assumption that Sbare{A,B,M]y) is local in the sense 
that an expansion in positive (but otherwise arbitrary) powers of derivatives exists, 
and that the terms of lowest order (without derivatives on B^^^, and F^iy) are non-zero. 
Its dependence on Myy is then dictated by power counting. 

Since the split of the path integral measures in ( |2.7|) -( p^ is somewhat formal, 
we will very briefiy present an alternative definition of the model. To this end we 



introduce in the path integral (|2.6|) a) sources J"^^^ and Jb^^u ^^ ^/^ ^^^ 5^^, and 
b) an infrared cutoff for A" and S" in the form of a quadratic term {~ASk{A, B)) 
in the exponent, with 

(2.10) 
The tensors R^ and i?"^ should suppress low momentum modes, i.e. become large 
or even diverge for p^ < k"^ (without having zero modes). First we choose k = M\y 
and define 



^-Wm^{Ja,Jb) 



[ VA VB VW VcDc (.-^YM-^SM^-Sgf-Sgh-SB+JA-A+JB-B 

Af J 

(2.11) 

The Legendre transform of Wmw^Jai Jb) with respect to J^, J^, the effective action 

rMiy(^,-B), plays the same role as Share{,A, B , Mw) in 



It is staightforward to replace My/ in ( p.ll|) by a varying IR cutoff A;, take the 



derivative of both sides of eq. (|2.11|) with respect to k, and derive an exact Wilsonian 



renormalization equation fl^ for Wu^Ja^ Jb) or its Legendre transform Vk{A^B). 
The A^ — B^y - model is then defined by the corresponding Wilsonian RGE for Wk 
or Ffc, and the boundary condition at fc = My/ given by ( p.ll| ). (The advantage 
of this procedure is that Slavnov- Taylor identities of the underlying SU{Nc) YM 



theory are well under control |T8| even in the presence of cutoffs as in (|2.1CI|) ). For 



convenience we will, however, most of the time refer to (|2.8|) as the definition of the 
model. 

Apart by locality, Sbare{A, B, Mw) in ( |2.8|) is constrained by the A^ = A^c — 1 U{1) 
gauge symmetries preserved by the second term in Sgf in ( |2.3| ), up to a covariant 
linear gauge fixing term corresponding to the first term in ( |2.3| ) (this can easily 
be proven by deriving a corresponding Ward identity). In the absence of fields 
carrying t/(l) charges Share can thus depend on A^ only through its abelian field 
strengths F^^,, hence we write Sbare{F, B, Mw)- However, Share is not constrained by 
renormalizability; we have to allow for arbitrary powers in the fields and momenta. 
Subsequently we assume that all dimensionful parameters are proportional to cor- 
responding powers of M^ (which also plays the role of an UV cutoff A = M^, cf 
and omit the corresponding powers of Mw 



Still, the most general dependence of Share on F and B does not yet allow to 
solve the model in a large A^ limit. To this end one has to assume that Share depends 
only on 0{N) invariant bilinear operators (quadratic in F oi B). This corresponds 
to the assumption that the low energy effective action of SU{N^ YM is dominated 
by quadratic SU{Nc) invariants; then 0{N), under which F" and 5" transform 
as A^-plets, appears as an accidental global symmetry. As we will show later the 
qualitative features of our results do, however, not depend on the large N limit and 
hence not on the global 0{N) symmetry. 

Let us define the following three Lorentz scalar bilinear operators 



N 



0,{x) 


n—1 


02{X) 


n — 1 


Osix) 


= 1:.b;^{x)b;^{x) . 



(2.12) 

n=l 

One of the simplest bare actions of the present class of models, which exhibits 
already all of the essential features, is given by 

ShareiF, B) = J d'x IChareiO,) + ^ {d.B;,)' + | {O.B;,)'] (2.13) 



with 

1 
2 
Chare in ( p.l3D may contain derivatives acting on the operators Oi defined in ( p.l2|) . 



B^iy = o ^t^'^P'^ Bp'^ ■ (2-14) 



Somewhat ad hoc we choose just bilinear "kinetic" terms for 5" ; generalizations 
are straightforward and will be discussed in section 5. An artefact of the ansatz 
( p.l3| ) with constant h, a will be the appearance of asymptotic states (poles in the 
propagators) with masses ~ M^/h; we do not expect the corresponding states to 
appear in an effective low energy theory of SU{Nc) YM, since these would be Ap- 
plets of 0{N). They could easily be avoided by replacing h, a in ( |2.13| ) by "higher 
derivative" terms 

/i(p2) , a(p2) with %2)^^(p2)^o for b'|>M^ (2.15) 



as is to be expected once these terms are generated by loops of the "charged" sector 
of SU{Nc) YM. Our subsequent main results will not depend on this choice, and 
unless stated otherwise we will continue to work with constant h, a for simplicity. 

In order to solve the model in the large N limit we have to make assumptions on 
the N dependence of the parameters in Share- These assumptions can be summarized 
by rewriting (|2.13| ) as 

Share{F,B) = jd'x \^N C .^re {^) + \ (S,^;,)' + | (9,i?;,)'| (2.16) 

where the coefficients of C'^are are independent of A^. 

Our subsequent aim is the computation of the (Euclidean) generating functional 
W{J) of connected Green functions, 



and the effective action obtained by the Legendre transform of W{J). The indices 
<Mw denote modes with momenta p^ < M^ as discussed below eq. (|2.7|). To this 
end we can add cutoff terms AS{A, B) to the exponent in ( 2.17|) , which have been 



defined in ( p.lO|) . Now, however, the cutoff functions R^ and R^ should suppress 
high momentum modes with p^ > M^. 

3 The large N solution 

The most convenient formalism for the solution of models of the form (|2.1(j|) in 



the large N limit is the introduction of auxiliary 0{N) singlet fields [jT9|, ^, which 
we will use as well. Already in the case of bosonic 0{N) models with spontaneously 
broken symmetry there are, however, some subtleties associated to the convexity 
of the effective action: In order to "see" the convexity of the effective potential in 
the broken phase it is necessary to introduce an artificial infrared cutoff k"^ (in, e.g., 
momentum space) and to study carefully the limit k"^ -^ 0. Conventionally this 
phenomenon is investigated in the context of the Wilsonian Exact Renormalization 
Group [^, but this formalism is not obligatory. 

In the case of models of the form ( |2.16| ) we observe a similar phenomenon: Now 
the "confining phase" is observed only in the limit of an infrared cutoff fc^ tending 
to zero, and again its existence is related to the convexity of the effective action. 



In order to review these points, which play an essential role in our class of models, 
we will first discuss them in the context of the bosonic 0{N) model in the broken 
phase in the next subsection 1. In the subsection 2 we will return to the A^ — i?^,^ 
model and use what we have learned before. 

3.1 The solution of the scalar 0(N)-model 

The field content of this model is given by an A^-plet of scalars ip"" , n = 1 . . . N. 
It is convenient to introduce a composite 0{N) singlet operator 

0{x) = V9"(a;)v9"(a;) . (3.1) 

The bare action of the model contains a bare potential, which can conveniently be 
expressed in terms of the operator 0{x). If one insists on renormalizability the bare 
potential contains only terms linear and quadratic in 0{x). Including its generic A^ 
dependence it is then of the form 

(We chose a minus sign in front of m^ in order to have m^ positive in the broken 
phase.) In order to implement both an infrared cutoff k and an UV cutoff A we 
supplement the kinetic terms with a cutoff term AS*^: 

SUv) = \l d^x 9^V^"9^v^" + /^St , (3.3) 

where 

^St = - j ^, V{p) R^,{p') V{-P) . (3.4) 

The cutoff function R^ diverges for p^ < k^ and p"^ > Is?, and vanishes for 
^^ ^ P^ ^ ^^- The aim is to compute the generating functional Wk{J) in the 
limit /c — > 0, where 



1 

~ 2 


//: 


X d^<p'' 


^9^V^" + /^st 


\l- 


d^P 

{2txY 


vij>) 


RUp') vi-p) 



The UV cutoff A is considered to be fixed, and A/"' ensures VFo(O) = 0. The auxiliary 
field method consists in writing 



10 






(3.6) 



In the large N limit the path integral in ( p.6| ) is dominated by the stationary point (s) 
and the relation between Vbare and Ghare becomes a Legendre transformation: 



NV^are ( ^ ) = NGl.arei'P) + j d^X <P ■ O 



(3.7) 



where = (piO) solves 



5(t){x 



-NGbare{(p) 



- 0{x) 



(3.8) 



HO) 



Conversely we have 



S fO 

60{x) '"'''[n 



${x) 



(3.9) 



Jew 



The correctness of the assignment of the powers of N in (|3.7| ) can be verified by 
rescaling O -^ NO', after which each term is proportional to A^ and the factor A^ 
can be dropped. 

Generally, for given Vbare, one can solve ( |3.9| ) for 0{(f)), insert it into ( p.7| ) and 
obtain GbareC4>)- With Vbare as in (|3l^) one obtains 



1 /" / 171 \ 

Gbare{4>) = ~T / '^^^ ( '^ + ^ 1 



(3.10) 



^ 77? A 



(3,11) 

Note that Gbarei4'i) is negative and unbounded from below for 0j -^ oo, which is no 
point of concern as long as is an auxiliary field with algebraic equations of motion. 
Using (p.6[) the expression (|3.5|) can be written as 



-WkiJ) 



- /p0 e"^^'""-'^^'^) f -p^ ^~s„Wp)~J dM<t>-o-j-v>n 



M 



(3.12) 



with M = M'M^. 

The Dip path integral is thus Gaussian and gives 



-Wk{J) ^ }_ I J)^ ^~NGtare{<l>)-NAG{<l>)+^ d*Xld^X2J''{xi)P{xi,X2,<l>)J"{x2) (S.IS) 



11 



with 



AG'(0) = --Trlog(P(xi,X2,0)) . (3.14) 

(A factor A^ from the trace in ( |3.14| ) has aheady been included in ( p.l3|) .) The 
propagator P{xi,X2,4>) has to be obtained from the quadratic terms in ip in the 
exponent in ( |3.12| ); it satisfies 



,4 S'^Skinj'f) 

Snk S'^iXi - X3) 



+ 2Snm S'^iXi - X2)(p{Xi) 



P"^^(X2,X3,0) 



(3.15) 



and we always use P™'^ = SmkP- 

For constant its Fourier transform reads 



p2 + R^(p2) + 2c/, 



and (p.l4|) becomes 



AG(0) = \jd'xj^^ log (p2 + Rt{p') + 20 



(3.16) 



(3.17) 



For what follows it is most convenient to employ a "sharp" cutoff function P^(p^), 
which is infinite for p^ < k"^ and p^ > A^ and zero for fc^ < p^ < A^. Up to irrelevant 
terms independent of this gives the same result as corresponding cutoffs of the 
d^p integral, and one obtains 



AG(0) 



d^x{ [a" - 402) log (a^ + 20) 

- [e - 402) log {e + 20) + 20 (a2 -e)} . 



647r2 



(3.18) 



Again, in the large A^ limit the path integral ( p.l3| ) can be replaced by the stationary 
point, and Wk{J) becomes 



Wk{J) = NGtarei^) + NAG{4>) 



rf^xirfVsl J"(xi)P(xi, X2, $)r{x2)} (3.19) 



where = 0(</) is the configuration which extremizes the right-hand side of ( p.l9|) . 
First we look for the vacuum configuration 0o = 0(J = 0) which has to solve, with 
Gbare from (|3A0|) and AG from (|318|) , 



12 



(G,a.e(0)+AG(0)) 



J 00 



->+'"''+IS5 



20 log 






+ K'-e 



(3.20) 



Clearly the solution 0o of ( p.20| ) depends on the infra-red cutoff k"^. It is instructive 
to follow the solution 00, as a function of the infrared cutoff /c^, from k"^ = A? down 
to k'^ -^ 0. For k"^ = A^, where AG(0) vanishes, we have 

2 

$0 = 3.21 

which is negative by assumption and our convention in ( p.2|) . For fc^ < A^ the 
solution 00 has always to correspond to a positive argument of the logarithm in 
( p.20|) , which implies (assuming A^ > m^, i.e. A^ + 20o > 0) 



Under the additional assumption 



^^>- 


2 


A2 
167r2 ^ 


m 
A 



(3.22) 



(3.23) 



the following subtle behaviour is observed in the limit of vanishing infrared cutoff 
k'^: One has simultaneously 



e^o , 






00 -^ 0_e , 






1 ^, fk' + 2$o\ 
Svr^'^^^^n A2 




A 


A2 
167r2 



finite 



(3.24) 



Note that the non-analytic behaviour ( p.24|) could not be obtained if one would solve 
( p.20|) right away for k"^ = 0. It corresponds in fact to the flat "inner region" of the 
effective potential Veff{<f"') in the broken phase, which is required by the convexity 
of the effective action (in terms of if"', once has been eliminated by its equation 
of motion) . 

In order to clarify this point we have to construct the /c-dependent effective action 
Tk{ip) from Wk{J) by a Legendre transformation: 



13 



Tkiv) = Wk{J) + / d^x r{x) (^"(x) 



f3.25) 



Taking the functional derivative of (|3.25| ) with respect to J"(x), and using the 
expression (|3.19|) for VFfc(J), one obtains for the "classical" field y?" 



^'\x) = -^^ = / d'x' P{x,x'A) r{x' 



(3.26) 



(j) in ( p.26|) is the stationary point configuration of Wk{J), as given by (|3.19| ), in the 
presence of sources: 



5(p{x) 



Scpix] 



iNG,are{<P)+NAG{<P)) 



+- / d'^xi d'^x'i d'^X2 d'^x'2 J^{xi) P{xi,x[, 

■P{X'„X2A) ^"(2:2) 



6P {Xi,X2,< 

6(j){x) 



(3.27) 



We have expressed the functional derivative of the term ~ JPJ in Wk{J) with 
respect to in terms of 6P~^/6(f) for two reasons: First, from ( p.l5| ), this functional 
derivative is trivial: 



Or [Xi,X2, 

(50(x) 



2S\x[ - x)5\x'2 - x) ; 



(3.28) 



Second, using ( |3.26| ) and ( |3.28| ), eq. (|3.27| ) becomes very simple once written in 
terms of (p"". one finds 
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6(f){x] 



(NGkarei^) + iVAGfe,,e(0)) + <^"(a;)^"(x) 



(3.29) 



Eq. (|3.29| ) allows to obtain as a functional of the classical fields (f"'{x) or 0{x) = 
ip"'{x)ip"'{x). From ( p.25| ) with (|3.19| ) for W{J), the inverse of (|3.26|) for J" in terms 
of if", and (|3.15|) for P^^ one obtains for the effective action Tki^p)'- 



Tkiv) = Sk^n{v) + N (Gbare{$) + AG{$)) + jd^X ^{x) ■ 0{x) (3.30) 



where cj) solves (|3.29|) . For the effective potential for fc = we replace cp"' by constants 
in (|3.30|) , i.e. we drop Skmif) and write 0{x) -^ O = y9"(^". Now we proceed to 



14 



solve (|3.29| ). For constant O the solution of (|3.29|) exhibits a similar behaviour in the 



limit of vanishing infrared cutoff as in ( |3.24|) : In the "inner region" of the effective 
potential, where 



^=lv^"r<^-^ (3-31) 



(instead of ( p.23|) ), one obtains 




^^-T7rT-^ = finite • (3-32) 

Hence, for \ip"\ satisfying ( p.31| ) we obtain = in the expression (|3.30|) for the 
effective potential, which is thus flat (yj-independent). 

As mentioned above, the underlying reason for this phenomenon is not the large 
A*" limit, but the convexity of the effective action of y?" even in the phase with 
spontaneously broken symmetry. However, its observation requires the introduction 
of an infrared cutoff k"^, and the study of the limits ( p.32| ). 



We discussed this phenomenon in quite some detail, since its essential features 
are shared by the A^ — B^^, - model, which is the subject of the next subsection. 

3.2 The solution of the A^ — B^^ — model 

Now we use a very similar formalism in order to solve the model defined by 
the partition function ( p.lTj ) in the large A^ limit. As discussed below ( |2.17|) , we 



may implement both UV cutoffs and infrared cutoffs A and k, respectively, with the 
help of cutoff terms ASk{A,B) as defined in ( |2.10| ) and appropriate properties of 
the cutoff functions R^ and R^ . (Again we keep the UV cutoff A ~ Mw fixed, but 
study the limit A; ^ 0.) Thus we write instead of ( [^.17D 



V' J 



WkiJ) _ 

(3.33) 



with Sbare{F,B) givcu by ( p.l6|) . 
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Next we introduce auxiliary fields 0j, one for each operator Oi appearing in the 
argument of C'bare in (|2.16| ) (recall the definitions (|2.12| ) for these bilinear Lorentz 
scalar operators). In analogy to eq. ( |3.6|) we write 



-Njd^xC'bare{^) ^ J_ f J)^. ^- NGtarei^^)- J d-^XlP.-O, _ (3.34) 



Again the path integral on the right-hand side of ( |3.34| ) can be replaced by its 



stationary point in the large N limit, and the relation between Ghare and C'bare 
becomes 

N J d^X C'bare (^) = iVGwe(<^^) + / d^x <j), ■ Oi . (3.35) 

In contrast to the scalar 0{N) model we do not assume here that C'bare contains 
only terms up to quadratic order in Oi. Hence the Legendre transform ( |3.35| ) may 
be quite involved in general, but this is not our point of concern: We may equally 
define our models right away through functionals Gbare{4>i)- Note that we allow for 
derivative terms in C'bare where the derivatives act on Oi] under the assumptions of 
locality specified below eq. ( |2.9| ) we will still obtain a local expression for Gbare{(pi) 
including, of course, derivatives acting on 0j. 

Replacing d^M ) into j^M) one finds (AT = AT' • Af^) 



(3.36) 

and the VADB path integrals become again Gaussian. The result can formally be 
written as follows: First we introduce the notation 

^",^: = {^;:,i?;.} , (3.37) 

i.e. the indices r, s, attached to the fields v?" denote both the different fields A", S" 
and the different Lorentz indices. Correspondingly we introduce sources J", J" for 
J^_^ and JB^fj.u- Now (|3.36|) becomes 
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with 



AG(0^) = -\ Tr log (P'-^(a;i,X2,0.)) 



(3.39) 



Again the propagators P^^ of the A" B^^-sysiem. are, in principle, proportional 
to 5n-m with n, 771 = 1 ... A^; we omitted these Kronecker symbols for simplicity and 
took care of the resulting contribution from the trace in ( p.39| ) by the explicit factor 
N multiplying AG in (|3.38|) . The propagators -P^,^, P^^pa and Pj^^p^ depend on 
the five terms A^^, 0^ ■ C^, ^{dBf, ^{dBf and \{dAf in the exponent of ( ^361) . 
Simple explicit expressions can be obtained only for constant fields 0j. They are 
given, in the Landau gauge a — > 0, in the appendix. For AG in (|3.39| ) one obtains 
for constant (/),■ 



AG((/) 
+4(01 



' 2J J (27r)4 
+ /?^(p2))(03 + i?^(p')) 



Iog(ap2(0i + i?^(p2)) 

+ log (/i/ + 4(03 + i?''(p')) 



(3.40) 



The cutoff functions R and R parametrize the simplest form of AS'fc given in ( [A.l 



in the appendix. Again we will employ, for simplicity, "sharp" cutoffs i?"^, R^ which 
are infinity for p^ < k"^, p^ > h? and zero for k"^ < p^ < h?. (We have checked that 
the stationary point in 0j described by eqs. ( ^.461 ) below appears also for smooth 
cutoff functions i?"^, R^ .) Then the ct^p integral in ( |3.40D can be performed with 
the result 



AG(0,) = 
+ (^A^- 



327r2 



d X 



, (40103-02)2^ 



A^- 



a 



2^2 



log (aV0i + 40103 - 



log(A;2(T0i +40 



1<P3 



(40103-02)2 - 



/ 0-01 



(A2 



log (A2/. + 403) - A:^ - ^ log {k^h + 40 



(3.41) 

This expression for AG{(j)i) has to be inserted into ( |3.38| ), and again the P0i path 
integral is dominated by its stationary point (s). Hence Wk{J) becomes 
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Wk{J) = NG{$,) -\j d^^i d^X2 {j;(a;i)P"(xi,X2>.)</r(a;2)} (3.42) 



where 



and (pi = (pii^J) satisfy the three equations (recall i = 1,2,3) 



(3.43) 



^{NG{4>i) -iJd'xi d^X2 J:ix,)P''{x,,X2,(Pi)J: 



.^2, 



. (3.44) 



■PriJ) 



First we are interested in the vacuum configurations 0° = 0i(J = 0). The three 
stationary point equations ( |3.44| ) with AG given by ( ^.41| ) simplify considerably if 
we switch from the three independent fields 



)1, 02, 03 to 01, 03 



S with 



nvz 



(J01 



(3.45) 



Then the stationary point equations for J = become (a factor N can be dropped) 



5G 



bare 



+ 



327r2 



(A- 



(3.46a) 



5G 



bare 



+ 



Att'^K^ 



ilog 



+ hk:^ 



403 + /iA2 



/i(a2-A;2)| 



, (3.46b) 



5G,are , 3 f^, ( ^ + k^' 



k^-k'' 



. 



(3.46c) 



Clearly eqs. ( |3.46D can not be solved without knowledge of Gbare (or its Legendre 
transform Chare, cf. ( p.35| )), and they may have several solutions. In the latter case 
the one with the lowest vacuum energy has to be chosen. This issue, together with 
the physical interpretation of the auxiliary fields 0j, will be discussed below. 

Let us now assume that, in some analogy to (|3.23| ) in the case of the scalar 0{N) 
model. 



5G. 



bare 



SE 



+ 



S0=0 



167r2 



A^ <0 



(3.47) 



Then, provided S°(/i;^ = A^) < 0, the solution of ( |3.46c| ) behaves as follows in the 
limit of vanishing infra-red cutoff: We have simultaneously 



(3.48a) 



t%e 



3 f.0, f^' + k'- 
L log 



IGtt 



SO + A^ 



6G 



bare 



5S 



S0=0 



167r2 



A^ 



(3.48b) 



(3.48c) 



From (|3.45|) the stationary point ( |3.48b|) , if it exists, corresponds to (assuming 0^ < 

oo) 



40^ 4>t~<P'2 =0 



(3.49) 



Eqs. ( |3.46a| ) and ( |3.46b| ) represent two additional constraints on the three fields 



(pi, 02 and 03, which are thus all determined. As in the scalar 0{N) model this 
stationary point would not have been observed if one puts k'^ = from the start. 

Subsequently we will denote the stationary point ( p3.48D , where ( |3.49| ) holds, the 
"confining phase" of the model. As stated above its existence depends on the validity 
of (|3.47|) ; clearly, this condition requires no "fine-tuning" of Gbare or Cbare- 

Next we construct the effective action of the model. Of course it is given by the 
Legendre transform of Wk{J) in (|3.42| ): 



UA, B) = w,{j) + J d'x {ji^a; + ji^^b;;) . (3.50) 

Following the steps discussed from eqs. ( |3.26|) to (|3.30|) in the context of the bosonic 
0{N) model one obtains 



T,{A,B) = AS, + NG{$,) + jd'x Uo. + | {d.B;,)' + | {d.B;^)') (3.51) 



with G(0i) as in ( p.43|) , and where it is again convenient to replace the sources in 
the stationary point equations ( p.44| ) by the fields A^, B^^ or the corresponding 



operators Oi ( |2.12|) : 
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-iVG(0,) + O, 



^<t>^{A,B) 



5V, 



(3.52) 



J MA,B) 



In the vacuum we require A" i?" = and hence Oj = in order not to break 



Lorentz invariance, and then p.52|) coincides with p.44|) for J = 0. 



Let us briefly discuss the relation between the condition TP{k'^ -^ h?) < for 
the confining phase to appear, and the convexity of the effective action in the case 
of a non-convex bare action. To this end we consider constant configurations F" , 
B^j^ in the expression (|3.51|) for Tk{A,B). 

Given the dependence of Oi on F^^,, i?^^, in ( |2.12|) it is straightforward to compute 
the determinant of second derivatives of Tk{Oi) with respect to F^^, B^^. If one 



drops ASk in (|3.51| ), this determinant is negative if 



401 03 - 0^ < . (3.53) 

In this case Tk without ASk would be non-convex around the origin F = B = in 
field space. Under the assumption 0i > (otherwise the term ~ F^ in Tk would 
have the "wrong" sign) ( |3.53| ) corresponds to the negativity of S, cf. eq. ( ^.451 ). The 
confining stationary point S -^ 0_£ for A;^ ^ corresponds thus to the "ironing" of 
a non-convex bare action as in the case of a non-convex bare potential in the bosonic 
0{N) model in the broken phase. 

The condition for the reach of the confining phase for A;^ ^ in the presence 
of no n- vanishing backgrounds A^, i?^^ and hence Oi differs somewhat from the 



positivity of the right-hand side of (|3.48c|) . After switching to the independent 
variables 0i, 03 and S, the three equations ( 3.52| ) can be written as 
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NG{<P,) + O 




(3.54a) 



-iVG(0,) + O3 




(3.54b) 



(5S 



NG{(Pi 




(3.54c) 



In principle we can solve ( |3.54a|) and (|3.54b| ) for 0i((9j,S) and 03((9i,S). Then 
equation ( p.54c|) has the "confining" solution S(A;^) -^ 0-e for A;^ — > if 



20 



Notably we then have, for A; — i> 



\ SE 167r2 J 202 

2 



< . (3.55) 

Js=o 



^103 - 02 



(3.56) 



even for non- vanishing classical fields A" i?" or Oi. Generally, once the right-hand 
side of ( |3.48c| ) is positive, there exists a finite region in the space of operators O^ 
where ( p.55| ) holds as well. 

Many remarkable properties of the effective action ( p.51| ) in the confining phase 



( p.56|) will be discussed in the next chapter. To close the present section, we will 
clarify the physical meaning of the auxiliary fields (pi which parametrize the vacua 
of the model. 

To this end we return to the partition function ( p.36| ). Instead of "switching on" 
sources J for the fields A"^, i?^^ we will consider sources Ki{x) for the composite 
operators Oi and define 

(3.57) 

(We could have kept the sources J; we have omitted them just for simplicity.) 
Evidently Wk{K) generates the Green functions of the composite operators Oi] 
notably we have 



One observes in (p.57 ) that the sources Ki couple to the bilinear operators Oi in 



the same way as the auxiliary fields 0,. Repeating the steps which have led to eq. 
( p.42|) for Wk{J) one thus obtains immediately 



Wu{K) = NG,are{<Pi) + iVAG (0, - K,) (3.59) 

with AG as in (|3.41| ), and where the fields 0^ = (pi{K) satisfy 
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Using ( p39D and (jpg) in (|38D one finds 



(3.60) 



J MK) 



<0,(x) >K = N 






(3.61) 



Hence in the classical limit (or for an infrared cutoff fc^ -^ A^) where AG vanishes, 
the vacuum expectation values of the composite operators Oi vanish as they should. 
In the limit of vanishing source Ki we can use ( |3.6CI| ) again to write 



< OAx) > 



K=0 



-N 



6G, 



bare\ 



(3.62) 



where the auxiliary fields (pi satisfy the stationary point equation (|3.60|) with K^ = 0. 
At first sight (|3.62| ) seems to be a trivial consequence of the Legendre transform 
( p. 351) , but now the (pi are the extrema of Gbare + AG including the quantum con- 
tribution. In any case ( |3.62|) shows that different vevs (pi parametrize different vevs 
< Oi > (recall their definitions ( p.l2|) ) . This relation would be linear if Gbare would 
be quadratic in (pi (as in the renormalizable 0{N) model). 

Since the (pi are auxiliary fields for composite operators the interpretation of the 
effective action (|3.51|) (at fc = 0, where ASk vanishes, and for constant field config- 
urations) as an energy density is not straightforward pO|, if the (pi are considered 
as independent variables. Notably T{(pi, A, B) is generically unbounded from below 
for \(pi\ — > cxD due to the negativity of Gbare{(pi)- (Such a feature of a potential V 
involving auxiliary fields is well-known from supersymmetric theories, once auxiliary 
F- or D-fields are introduced in order to complete the supermultiplets. Then one 
has e.g., V{F,ip) = — l-Fjp + FiWi + h.c. where Wi is the derivative of the super- 
potential W{ipi)). At the stationary points (pi = (pi, however, the interpretation 
of T[(pi, A = B = 0) as an energy density can be maintained [^ and the energy 
densities of different vacua can be compared. 

The interesting properties of the effective action in the confining phase with 
( p.49|) will be the subject of the next section. 
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4 Properties of the confining phase 



In the following we assume that the necessary condition ( |3.55|) for the existence 
of a confining phase is satisfied within a finite regime of field configurations A^, 5" 
(and corresponding bilinears Oi), and that it represents the lowest (or only) vacuum 
of the model. For k'^ —>■ the effective action ( p.51D reads 

T{A, B) = NG{^.) + / d'x ho. + ^ {d,B;:f + | (a,S;,)'| (4.1) 



with G{(j)i) as in ( p.43|) , and the 0j satisfy (|3.52| ). The confining phase is characterized 
by the validity of eq. (|336|) , 40i03 -$l = 0. 

First we consider a new kind of N U{1) gauge transformations which involve 
vector-like gauge parameters A^ 0: 



6A;ix) 



5b;{x) 



a;{x), 5f;^{x) = d,K{x) - dAl{x) 



A"(^) 



\ 



I KA-) 



(4.2) 



Using eq. ( p.56| ) in order to eliminate 02, (choosing 02 positive, otherwise the sign 
in 5B has to be changed) and ( p. 121 ) for the definitions of the bilinear Oi, one finds 
that the terms (piOi in (14.11) can be written as 



^iO, = E 



H f;^ + 



^3 5 



HV 



(4.3) 



Hence these terms are invariant under the gauge transformations ([4.2| ). Next we 
have a look at the term ~ (c)^S" )^. Thanks to a Bianchi identity one obtains 



^id,B;^ 



+ o a 



'/iVi 



(4.4) 



i.e. it is also invariant up to derivatives acting on 0j which appear in the variation 
55" in ( [4.2|) . Recall that the (pi are fixed in terms of the bilinear operators Oi 



by (|3.52| ); for constant configurations Oi the 0j would also be constant, and the 



symmetry breaking contribution to ( ^4.4| ) is thus proportional to d^Oi. 

The last term ~ {d^B^^Y in (^?T| ) behaves as a 
([4.2|), and altogether one obtains a Ward identity 



The last term ~ [d^B'^^Y i^ (jl.lp behaves as a gauge fixing term of the symmetry 
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6T{A,B) = -, 



\ 03 



(4.5) 



It is natural to expect that the syinmetry ( [4 .21 ) is broken by higher derivative 



terms (beyond the gauge fixing term): The bare action Share ( P-16| ) of the model 
does certainly not exhibit the symmetry ( ^.21 ) (or satisfy the Ward identity ( ^4.5| )), 
and the Green functions at large non-exceptional Euclidean momenta with p^ — > A^ 
are generated by Share- This fact is realized by the dependence of the effective 
action on higher derivative terms. The symmetry ( [4. 2] ) is thus a pure "low energy" 
phenomenon. 

Note, however, that the "violation" of the Ward identity ( [4.5|) is of higher order 
in the fields, and is proportional to derivatives acting on the bilinears Oi. Hence the 
symmetry ( [4.2|) is realized on modes A^ and 5" which may be rapidly oscillating, 
but correspond to constant bilinear configurations Oi. (These modes are associated 
with Green functions with "exceptional" external momenta.) 



The implication of the gauge symmetry ( [4.2| ) on these modes of the U{1) gauge 
fields A^ is evidently that they can be "gauged away" and "eaten" by the (massive 
or even infinitely massive) Kalb-Ramond fields -B" 0. 

Correspondingly the (approximate) gauge symmetry ( |4.2| ) will have consequences 
on the "response" of the model with respect to couplings to external sources J^ 
and J^ . Before studying this issue we will discuss the behaviour of the effective 



action ( |4.1| ) under a duality transformation, whose existence is closely related to the 
symmetry ( |4.2| ). 

Let us recall the essential features of a duality transformation: It corresponds to 
the introduction of dual fields such that the Bianchi identities associated to the dual 
fields correspond to the equations of motion of the "original" fields (up to gauge 
fixing terms), and the Bianchi identities associated to the "original" fields to the 
equations of motion of the dual fields. 



Using (|4.3|) , the equations of motion for A^' and B^,^ are 



d, [<PiF;, + \J<PihB;,] = , (4.6a) 

2^,F;^ + 203i?;. - ^e,.,.dpd^Bl^ - I (d.d.B:^, - d^d^B:^) = . (4.6b) 
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As dual fields we introduce vectors C!^ (whose field strength tensor will be de- 
noted by F^i^), and pseudoscalars </?". They are related to the original fields A^, 
B^j^ through 



_ pc,n 

2 f"^ 



/lU 



^iK. + ^MzBl 



fll/ 



(4.7a) 



d.v'' + c: 



h 



2\ 



d,B;^ 



(recall our definition of B'^^^ etc. in eq. ( p.l4| )). 

The Bianchi identities associated to the dual fields are 



(4.7b) 



d,F%^ = 



(4.8a) 



£^,upa dpd^if'^ = 



(4.8b) 



Contracting the dual of ( [4.7aD with 9^ one realizes immediately that the Bianchi 
identity ( [4.8a| ) coincides with the equation of motion ( |4.6aD . The Bianchi identity 
( [4.8b| ) should reproduce the equation of motion ( ^.6b| ) up to the "gauge fixing term" 
~ a. Indeed this is the case up to terms proportional to (9^0i. 

The Bianchi identities associated with the original fields A" 5" read 



^mC = 



(4.9a) 



d^d,B;^, = . (4.9b) 

As action for the dual fields C" cp^ we propose (up to gauge fixing terms) 



Dual 



I 



d X 



nv3 



/lU flU 



+ 



I h(B uj^' + r 



(4.10) 



Now one finds, using ( [4.7b| ), that the Bianchi identity ( |4.9bD coincides exactly with 
the equation of motion for yj" derived from ( [4. 101) , whereas the Bianchi identity 
( [4.9aj) agrees with the equation of motion of C^ only up to terms ~ 5^0j. 
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Not astonishingly, duality is thus realized to the same level as the gauge symme- 
try (|4.2| ): To the order quadratic in the fields the physical content of the dual action 
( [4.10|) is the same as the original effective action ([4.1|) , since the "perturbing" terms 



~ d^(j)i are effectively of higher order. (As solutions of the stationary point equations 
( p.52|) the deviations of the fields 0j from their constant values in the vacuum are 
proportional to the bilinear operators Oi.) To higher order in the fields duality is 
realized only on modes of A" i?" which correspond to constant bilinears Oi. 

The physical interpretation of the dual action ( [4 .101 ) is obviously the one of 



an abelian U{1)^ Higgs model in the spontaneously broken phase where ip"^ rep- 
resent the Goldstone bosons, and where the gauge fields C^ have acquired a mass 
2(03//i)^/^. Since this represents the "low energy effective action" of a theory in 
which the "dual" electric charge has condensed in the vacuum, the original action 
( |4.1D corresponds to the situation where the "magnetic" charge has condensed in 
the vacuum. 

Let us return to the response of our model with respect to external sources. The 
expression for W{J) has been given in eq. ( p.42|) in the preceeding section, and 



subsequently we consider the limit of vanishing infrared cutoff fc^ — > 0. We recall 



that the fields 0j(J) satisfy the stationary point equations ( |3.44| ), and satisfy thus 
obviously 

UJ) = ^' + 0{J') (4.11) 

(where the 0° satisfy (|3.49| )). 

In fact, the procedure to solve (|3.44|) exactly would be quite cumbersome, since 
one would need the expressions for the propagators P'''*(xi,a;2, 0i) for arbitrary (pi. 
In practice it is much easier to work with the effective action T{A.,B) (|3.51| ); this 
procedure will be described below. 

On the other hand, if the sources J are sufficiently weak, we can expand W{J) 
in powers of J . Since the 0° are stationary points of G{(f)i) we have 

W{J) = NG{^^) - \ fd^xi <fx2 {Jr(^i) P'-^(xi,a;2,0°)Jr(x2)} + 0{J^) . (4.12) 



Interesting informations can already be obtained from the second term of 0{J'^) 
in eq. ( [4.12|) , which involves the propagators P^^ in a confining vacuum given in 
eqs. (|A.7|) and ( |A.8|) in the appendix. 
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Let us start with a source Jaii{^) fo^^ the fields A" only. The simplest geomet- 
rical configuration is a Wilson loop source J^ which is non- vanishing only on a 
curve C embedded in 4d space-time, where C has to be closed because of current 
conservation: 



Jl^{x)=igA dx'^6\x-x') . (4.13) 



The term quadratic in J in ( [4.12| ) then becomes 

with P^^ as in eq. ( [A.?] ). The divergent constants in ( |A.7| ) disappear actually in 



the expression ( |4.14| ). This can be seen most easily by applying the (abelian) Gauss 



law, and expressing the source J\ in terms of a source J"^ ^^^ which is non- vanishing 
on an (arbitrary) surface S bounded by the contour C: 

Jl^{x) = 2d.Jl,^{x) (4.15) 

with 



Jl,u{^) = Wf / d'a^^iz) 6\x - z) (4.16) 

where gp = Qa/'^i z parametrizes the surface S and 

d^a^y{z) = e^upa nl{z) n\{z) d^z . (4.17) 

n^p, nl are two orthogonal unit vectors perpendicular to S. 



The expression ( [4.14|) then turns into 



2Ngl f d'a,.{z) f rfV,.(^') P/^J^^^iz - z') (4.18) 



with P^^ as in eq. (|A.8a|) . The dependence on the divergent constants in eq. (|A.7|) 



has disappeared as announced. 

In the limit where the surface S becomes very large one finds that the expression 
( [4.18|) is proportional to the (minimal) surface bounded by the contour C which we 



will also denote by S. 

Thus we have obtained the area law for the expectation value of the Wilson loop. 
This result has to be swallowed, however, with two grains of salt. First, in the limit 



where S and hence the expression ( [4.18| ) become large, it becomes inconsistent to 
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confine oneself to terms of 0{J'^) in W{J). (Fortunately, as we will discuss below, a 
more complete treatment just modifies the contributions to the string tension, not 
the area law.) 

Second, and more importantly, an inconsistency arises if we couple our model to a 
quantum field theory. In order to clarify this inconsistency we first consider a slightly 
modified configuration of sources Jp'- Let us consider two distinct "Wilson surfaces" 
5*1, 5*2, which may both be small and centered at Xi, X2, respectively. Accordingly 
the source Jp becomes the sum of two terms, Jp = Jp^i + Jp^2, where Jp^i and Jp^2 
are non- vanishing only near Xi or X2- Inserting Jp into the term quadratic in J in 
W{J) in ( |4.12| ) one obtains, apart from "self-contractions", a mixed term which is 



of the order of -P^^o-(^i ~" ^2) for \xi — X2I large compared to the diameters of the 
Wilson surfaces. (The Lorentz indices depend on the orientations of the surfaces 5*1 
and S2.) Due to the second term in ( |A.8a|) P^^{xi —X2) decreases only as \xi—X2\~'^ 
for generic relative orientations. 

The resulting expression can be interpreted as a contribution to the action of 
the configuration: The effective action F is given by the Legendre transform ( p.50|) . 



and up to 0{J^) the expression for F thus coincides with the expression ( [4.12|) for 
PF up to a change of the sign of the second term. 

At this point the result is not yet problematic. Let us now consider the cou- 
pling of our model to a quantum field theory in the form of Jp uyFuvi where the 
source term Jp^^^ is quadratic in "external" quantum fields. {SU{Nc) Yang-Mills 
theory, e.g., contains couplings of the form fnabFJ];^W^Wl^, which corresponds to 
Jf,ij.u = fnabWI^Wl^.) Vacuum bubbles of the quantum fields then correspond to a 
"background" of sources Jpn^, which fills the whole space-time (to be averaged over 
•^F,fiu with, e.g., a Gaussian measure). If the contribution to the action due to the 
induced interactions between two vacuum bubbles decreases only as the (distance)"^ 
as described above, space-time filled with vacuum bubbles will lead to an infinite 
action. Even more severely, every additional Wilson loop switched on "by hand" 
will also lead to an additional divergent contribution to the action due to its inter- 
actions with the vacuum bubbles (after averaging over Jp^,, this divergence is only 
logarithmic at large distances). Hence it seems to be disallowed^. 

^Notc the important difference to ordinary QED: Here tlie contribution to the Euciidean action 
induced by two vacuum bubbles, i.e. two dipoles locahzed both in space and time, decreases as the 
(distance)"^. "Switching on" an additional single dipole immersed in a bath of disoriented dipoles 
then only costs a finite amount of action. 
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Actually, before we have assumed that we have obtained our model as an effective 
low energy theory after the fields W"" have been integrated out. The present paradox 
persists, however, if previously we have only integrated over PF^-modes with p^ > k^ 
(with k^ as small as one likes, but 7^ 0), and the remaining modes with p^< fc^ 
generate the vacuum bubbles described above. In addition, in full QCD we have 
more fields like quarks, whose vacuum bubbles will generate a background density 
of sources Jp^^^- 

The fact that we assumed that PV^-modes with p^ > A;^ have already been 
integrated out will, on the other hand, lead to a solution of the problem: First, 
in the original SU{Nc) Yang-Mills theory, a gauge-invariant Wilson loop "source" 
is coupled both to A^ and to VF^. If we integrate over W^-modes after having 
introduced the fields 5^^, as in eqs. ( p. 51) and ( p.6|) , this will necessarily generate a 
coupling of the Wilson loop source to B'^^. Hence it is inconsistent to switch on a 
source for F^^^, but not simultaneously a source for 5" . The simplest expression for 
a source for i?" is again of the form (|4.16| ) 

JbM^) = WB f d'a,,{z) 6\x - z) , (4.19) 

where the surface S coincides with the one in ( |4.16| ). The value of the coupling 
constant qb can, a priori, not be predicted within our model (neither the couplings 
gp or QA in (IH) or (glSD). 

Thus we return to the term quadratic in J in (|4.12| ), and insert for the source 
J a sum Jp + Jb with Jp as in ( [4.16| ), and Jb as in ( [4.19|) . With the propagators 
from (|A.7|) and (|A.8|) this term becomes 



Ibn'^a Js Js ^ ^ V 0i 

+0{\z,-Z2\-')} . (4.20) 

Thus the terms of 0{\zi — ^2!"^) in the curled parenthesis in ( [4.20|) vanish if 



^3 Qf = \l<Pi 9b ■ (4.21) 

Again, at this point this is not an obligatory constraint. However, if we switch 
on a) 'Wilson loop" sources Jp, Jb, b) simultaneously uncorrelated "background 
sources" Jp, Jb, c) average over Jp, Jb with an arbitrary Gaussian measure, we 
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obtain a logarithmically infrared divergent expression, unless the condition (^4.21 



on Jp, Jb is satisfied. Hence, in the background of a fiuctuating vacuum it costs 
infinite action to turn on sources Jp-, Jb which are not related through ( [4.21|) . This 



constraint can be interpreted as a constraint on the ratio (f)i/4>3 in the presence 
of sources, or on the ratio Qf/qb (once gp and Qb vary, e.g., with the vacuum 
configuration of a "microscopic" theory) or, most realistically, on a combination of 
both. 

Before discussing the resulting behaviour of the expectation value of the Wil- 
son loop we note that the above consideration leads, in general, to the following 
constraint on sources Ja, Jb or Jp, Jb'- 

JlAx) = 2^, d^Jlix) or (4.22a) 



^3 -JaA^) - ^\JVl '~'y^B,ufi 

'^^ JfA^) = ^i Jl,Ax) . (4.22b) 

The constraints (|4.22|) are obviously related to the gauge symmetry (|4.2| ): If we 

corre- 



replace (pi by the vacuum configurations (in the absence of sources) 0j, (|4.22 



sponds to the "current conservation condition" , i.e. to the condition that 

Jd'x{jl^A; + Jl^,B;J) (4.23) 



is invariant under ( [4.2| ). 

In the case of "conventional" gauge symmetries the current conservation condi- 
tions can (and have to) be imposed by hand in order to ensure renormalizability and 
unitarity of the theory. In the present model, on the one hand, current conservation 
cannot be imposed from the beginning, since the associated (approximate) gauge 
symmetry appears only at the level of the effective action once the equations of 
motion of the fields (pi are satisfied. On the other hand, renormalizability is not an 
issue here, since we consider an effective low energy theory with a fixed UV cutoff. 
Also unitarity is trivial as long as we assume the absence of poles in the propagators 
in the Minkowskian regime, cf. the discussion before eq. (|2.15|) . Thus the model 
has no S'-matrix at all, i.e. no asymptotic states. (The absence of bound states will 
briefly be discussed below). Consequently at this level no constraints on the sources 
arise. 

Only if we couple the model to a quantum fleld theory we have to reconsider 
the question of unitarity, i.e. the possibility to project - in a Lorentz covariant 
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way - onto a positive semi-definite part of its Hilbert space wliicli is represented, 
loosely speaking, by our sources J viewn as functionals of fields. Precisely in this 
situation "vacuum bubbles" impose eqs. ([4.22|) , the analog of current conservation 
in an abelian Higgs model with an additional Lorentz index attached to the currents. 



If eqs. ([4.22|) hold the expectation value of the Wilson loop has to be reconsid- 



ered: The Wilson loop corresponds now to a source J a of the form ( |4.13|) , which 
can be rewritten as a source Jp of the form ( [4.16| ) using ( f4.15| ), plus a source Jb of 
the form ( [4.19| ) where the surface S coincides with the one in ( ^.16| ). 

First, the term quadratic in J in ( [4.12[ ) consists now only of the "short range" 
contributions neglected in ( |4.2CI| ) . Using again the propagators ( [A.7| ) and ( |A.8| ) these 
read 



/ d^a^^^izi) I d^ap„{z2) 




IGvr^^i V h 



X Tr,,u,p.{d) 



h 



>Tc 



2,^iu,pa 



\Z\ - Z2\ 



r^i 2|zi-2;2 




(4.24) 



With ( [4. 211 ) the contributions ~ \zi — Z2\~'^ cancel as well, and we recall that the 
tensors Ti and T2 are given in ( [A.2| ) in the appendix. 

In the limit where the surface 5* becomes large the expression (^.241) behaves as 




dz Ki 




(4.25) 



Hence it implies the area law in spite of the cancellations of the long range contribu- 
tions of the propagators. (Since we have omitted the UV cutoffs in the space-time 
propagators the expression ( |4.25|) is seemingly UV divergent.) 

Expressions of the form ( [4.25| ) for the (negative) exponent of the expectation 
value of the Wilson loop have already appeared repeatedly in the literature in the 
context of the method of vacuum correlators 0, ^: In the Gaussian approximation 
(the stochastic vacuum model) the expectation value of the Wilson loop is given 
by the expectation value of the field strength correlator, and many models for this 
correlator lead to ( [4.25] ) |^, |T^ . Actually already the previous result ( [4. 181 ) - the area 
law for Wilson loop sources Ja only - has an interpretation in this approach: It would 
correspond to a function Di(x^) in the standard decomposition of the field strength 
correlator f^, ^ ^ which decreases only as \x\~^ for large |a;|; such a behaviour also 
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implies the area law, but it is strongly disfavoured both phenomenologically and 
from lattice data |2^ . 



As we have already stated several times above, sources for operators in a SU{Nc) 
Yang- Mills theory correspond in our model - if considered as an "effective low energy 
theory" - to a priori unknown superpositions of sources for A'^ (or FJ^^) and -B^^,. 
The previous discussion leading to eqs. ([4.22| ) fixes this ambiguity, and precisely 



in this case the two-point correlators decrease exponentially in agreement with the 
SU{Nc) lattice data. 



In our model the expression ( |4.25| ) for the resulting string tension can, however, 
not be taken too seriously, since in the limit of a large surface S the terms of higher 
order in J in W{ J) in ( ^.121 ) can no longer be neglected. 

In order to compute W{J) in the "non-linear" regime one has to start first with 
the effective action T{A, B) as given in eq. ( p.51|) (in the limit of vanishing infrared 
cutoff k) or eq. (|4.1| ). One has to solve the combined equations of motion for A" 
i?^j, and (pi in the presence of sources J^^^ and Jb^^u'- The corresponding equations 
for A" and 5" are given by eqs. ([4.6|) (where, however, the zeros on the right-hand 



sides have to be replaced by the corresponding source terms and 2y0i03 has to be 
replaced by 02)- The equations for 0j are given by eqs. ( p.52| ) or ( |3.54| ), where the 
operators Oi depend on F^^, and -B^j, as in ( |2.12| ). 



The corresponding fields have to be inserted into T(A, B) in ( ^4.1| ), and then one 
has to "undo" the Legendre transformation ( |3.48| ) in order to obtain W{J). As 
already noted this last step is actually trivial since T[A, B) is quadratic in A, B: It 
suffices to change the sign of the second term in the expression ( |4.1| ) . Then one can 
study the dependence of W{J) on J in its full beauty. 

Albeit generally the solution of the combined equations of motion is certainly 
quite involved, there is one situation where it becomes straightforward: Let us 
assume that we have turned on sources Jp^^ui^) and JB,^ui^) which are related as 
in eq. ([4.22b| ), and which are constant inside a space-time volume V, and vanishing 
outside. Let us furthermore assume that we can neglect all derivative terms in 
G{(f)i). Then we can allow for constant configurations (fi[ inside the volume V , and 
constant (generally different) configurations 0° outside V] the discontinuities at the 
boundaries of V cost no energy in this case. The arguments leading to eqs. ([4.22|) 



are based on the long distance behaviour of the propagators in the vacuum outside 
V , accordingly eqs. ( |4.22| ) have to hold for 0° outside V . For simplicity we assume 
that they also hold for (j)\ inside V , i.e. that (pX /(j)\ = 0?/03- 
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For the fields FJ^j^i^x) and B^j^{x) the following simple solution is then possible: 
For BJ^^{x) one chooses 

B;^{x) = (4.26) 

everywhere, both inside and outside V (because of the derivatives acting on B^^^ in 
the action discontinuities in B are not allowed). Then eqs. ( [4.6|) - with source terms 
included and using the "confining" relation 02 = 2-^/0103 both inside V and outside 
V - collapse to 



2dM,F::^{x)] = d,J^^^,{x) , (4.27a) 



khF^^x) = Jl^^ix) . (4.27b) 



With Jp and Jb are related as in (|4.22b|) inside V eqs. (|4.27|) are trivially solved 
simultaneously by 

C(^) = 4 JF,,uix) inside V , (4.28a) 

01 

FJ^^{x) = outside V (where Jp = Jb = 0) . (4.28b) 

Hence, in this simple scenario the effective action ( [4.1|) is given by 

r= f d*x {nG{$Y) + $iOi} + f d^x NG{$^) (4.29) 

where G^fpi) is the density associated to G{(j)i), and V denotes the entire space- 
time volume outside V. (Recall that only Oi is non-zero for 5" = 0.) Since the 
configurations 0° are local maxima of G{(j)i) and the configurations 0]^ are local 
maxima of G{(j)i) + 0iOi one derives easily 

iVG(0r) + 0rOi > iVG(0°) . (4.30) 

Hence the presence of constant sources inside a volume V increases the Eu- 
clidean action by an amount which is proportional to V (in an approximation where 
derivative terms in G{(f)i) are neglected). 
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This consideration can be applied to the Wilson loop - i.e. sources Jf, Jb given 
by ( [4.16|) and (|4.19| ) with gp and qb related by y^^gp = y(t>i 9b ^ provided we 



replace the enclosed surface S* by a "layer" of finite thickness: In the presence of an 
UV cutoff A we cannot resolve distances smaller than A~^ anyhow, and hence we 
should replace the surface 5* in ( [4.16| ) and ( [4.19[ ) by a layer of thickness A^^ and a 



corresponding volume V = S ■ A~^ . 

Given the above discussion one then obtains immediately an action F which is 
proportional to S, and - after a Legendre transform - W{J) proportional to S and 
thus the area law. In the present scenario the profile of the flux tube associated to 
the Wilson loop would be discontinuous, with 0j = (pf (and F^^, given by ( [4.28a| )) 
inside the flux tube, and (pi = 0° outside. If we would take derivative terms in G(0i) 
into account, this profile would be smoothened, and additional contributions to the 
action proportional to the length of the loop C would be obtained. As stated above, 
a treatment of these effects would require the solution of the combined equations 
of motion of (pi, F" and i?" (in some similarity to the much simpler scenario 
based on the dual abelian Higgs model) and depend evidently on the unknown bare 
action Share in (|3.33|) or Gbare in ( |3.35|) , ( p.37|) . 



In order to close this section we briefly report on our fruitless search for bound 
states in the model in the conflning phase. In principle these would show up as poles 
in the propagators of the flelds (pi for g^ < 0. To this end one has to develop the 
fields (pi around the vacuum configuration, 



(P^ix) = (Pi + 5(Pi{x) . (4.31) 

Inserting (|4.31| ) into the action ( |4.1| ) the propagators for (pi are obtained from the 



terms quadratic in 5 (pi. It is sensible to study first the "mass matrix" Mf, in the 



space 6 (pi: 

N 5''G{(P;) 



V 2 



(4.32) 
We recall that (7(0^) is composed out of G^are and AG, cf. eq. (|3.43|) , with AG given 



in ( ^.41 ) for A;^ — > 0, and the 0° satisfy ( |3.49| ). Due to the non-analytic behaviour 



of AG in the confining phase one finds that all entries of Mf, actually diverge for 
A;^ ^ 0, but Mj^- has two finite eigenvalues (depending on Gbare) with eigenvectors 
given by 
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V 



5? \ 



/ 



( \ 



vv^w 



(4.33) 



We have computed the full momentum dependence of the contribution (AP^^(g^))a5 
to the inverse propagators, which arises from the quantum contribution AG ( |3.39|) 
after expanding (pi to second order in 50j in the "directions" specified by ( |4.33|) , 
for \q^\ < A^. We found that {/S.P^^{q^))ab depend only weakly (logarithmically) 
on g^, and that its determinant is negative definite. Together with the negativity 
of Gbare{<Pi) (rccall the remark below eq. ( 3.35 )) this implies that the full inverse 



propagators do not vanish for \q^\ < A^, thus no bound-state pole appears at least 
below the UV cutoff. 

Such a pole would actually be a disaster for unitarity: The negativity of G{(f)i), 
notably for |0j| -^ oo, is harmless only if the fields (pi are pure "auxiliary fields", i.e. 
with algebraic equations of motion. 

On the other hand one may be deceived by the absence of "glueballs", if the 
present model is interpreted as an effective low energy theory for SU{Nc) Yang- 
Mills theory. However, we recall that the present model would only describe the 
abelian subsector of SU{Nc) Yang-Mills in the MAG, and that we finally have to 
add the (massive) "non-abelian" gluons W^ as well. Our theory induces flux-tube 
like confining attractive interactions among all fields which couple like JfuuJbuu 
to F and 5, hence the VF^-gluons will necessarily form bound states which will 
correspond to the desired glueballs. (This scenario is quite consistent with heavy 
glueballs, given Mw ~ C(l GeV) in the MAG 11). 



Further details on the possible relation between the present model and SU{Nc) 
Yang Mills will be discussed in the next section. 

5 Discussion and Outlook 

In this final section we want to discuss some general features of the model, 
notably the relation between the large N limit and SU{Nc) Yang Mills at large Nc, 
and the systematic inclusion of higher derivatives in the "bare" action. 



First, we had emphasized the solvability of the A, 



B^y - model for A^ 



cxD m 



the sense that W{J) or T{A^B) can be given explicitely for arbitrary Gi,are{4>i) as 
in ( p.42| ) or ([4.1|) , and one is left with the coupled equations of motion for A" 5" 
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and (pi. To this end we had to assume a certain A^-dependence of the coefficients of 
the bare action Sbare{F,B), cf. eq. ( p.l6|) : 



Share{F, B) = N Share ( ~7m' ~7m I ^^'^^ 



where the coefficients of Share are A^- independent. 

Let us now assume that S},are{F,B) has been obtained from SU{N^ Yang-Mills 
after integrating out the "non-abelian" gluons W^ as in eq. (|2.9|) . From a sim- 
ple analysis of the contributing Feynman diagrams and employing a^ ~ A*""^ one 
obtains, in the large N^. limit, 

SlfeiF, B) = Nl Slti f 4^, -4r] (5.2) 



in contrast to ( |5.1|) . In principle it is still possible to introduce the operators Oi 
( p.l2|) , and to express exp{—S^arei^i)) i^ terms of a path integral over (pi as in 



(PI: 



-S,^ii(0,) ^ 1_ fj)^. ^-N^Gi,rei^,/N,)-Jd^x^,0, _ ^5 3^) 



The powers of N^ on the right-hand side of ( [5. 3D have been introduced such that, 
in the stationary point approximation, the coefficients of Gbare are of 0{1). (This 
becomes clear after a rescaling (pi — > Nc(pi, Oi — > NcOi and using (|5^ )). However, 
now the stationary point approximation of the path integral (|5.3|) can no longer be 
justified in the large Nc limit. 

In addition, the structure of the operators Oi (|2.12|) assumes a global 0{N) sym- 



metry which corresponds to a rotation in the space of the A^ = Nc — 1 U{1) subgroups 
of SU{Nc) (or among the fields A" -B" , n = 1 . . . N). This 0{N) symmetry is not 
a subgroup of SU{Nc) and is generally not a symmetry of S]^Jf^{F, B). (It appears 
as an accidental global symmetry, however, if one assumes that S^Jf^{F, B) contains 
just quadratic SU{N^ invariants.) In general it is thus necessary to introduce A^ 
copies of each operator Oi in ( p.l2| ), 

Otl=F;,F;^ etc. forz = 2,3 (5.4) 

without sum over n on the right-hand side (or A^ independent linear combinations of 



36 



Cj- ) I . In particular for small N the introduction of operators Oi" is straightfor- 
ward, but again the stationary point approximation of the path integral operators 
Oi ( |5.3|) can no longer be supported by a large Nc limit. 

We emphasize, however, that the essential phenomena described by the present 
class of models do not depend on the large A^ limit (as the phenomenon of sponta- 
neous symmetry breaking in bosonic 0{N) models): These are the appearance of 
a confining phase, associated with an approximate (low energy) vector gauge sym- 
metry and duality transformation, and associated with the area law for the Wilson 
loop. They have their origin in the convexity of a quantum effective action in spite 
of a non-convex bare action, which holds independently of a large N limit. Just the 
precise relation between the effective action and the bare action can only be obtained 
for A^ -^ oo. An application of the present approach to SU{Nc) Yang Mills will thus 
require - apart from a determination of 5'^re ^s in eq. ( |2.9| ) - a performance of the 
X>0j path integrals beyond the stationary point approximation. 

In general we cannot expect that a "bare" action as obtained by eq. ( p.9|) depends 
only on operators Oi of the form (|2.12|) , since it will involve the most general terms 
with higher derivatives. On the one hand we recall that the functional Gbare{4'i) 
introduced in (|3.34|) , which specifies the models solved in Section 3, may well contain 
arbitrary derivatives acting on 0j. These functionals Gbare{4>i) correspond to bare 
actions Sbare{Oi) with arbitrary derivatives acting on Oi. On the other hand, these 
do not correspond to the most general form of higher derivative terms: In general 
it will be necessary to consider bilinear operators of the form 

(and with F replaced by B) where / depends on Kronecker symbols or ^-symbols 
in the Lorentz indices and derivatives d^ (with Fid^F2 = Fid^F2 — {d^Fi)F2). 

Nevertheless it is possible to proceed as before: One can introduce auxiliary 
fields 0i.Ai.-Ai for each of the operators in ( |5.5|) , and rewrite the most general bare 
action Sbare{Oi) as 



e 



-StareiOi) _ _ [ J) A . ^-Gtarei^i)- J d'^X(j>iOi /g Q\ 



■^Note that a Z2 symmetry (charge conjugation of all W^ together with A^ — > —A^) is a 
subgroup of SU{Nc), which allows to restrict ourselves to bilinear operators. 



37 



(adding, if one wishes, appropriate factors of A^ as in ( |3.34|) ). Of course the appro- 
priate contraction of Lorentz indices is understood in (piOi. 

Inserting (|5.6|) into (|3.33|) one obtains an expression for Wk{J) analogous to 
( p.36| ), and the VADB path integrals are still Gaussian. The essential technical 
complication is now the construction of the propagators P^** of the A^, B^^, - system 
which appear in (|3.38|) , ( |3.39|) , and the computation of AG'(0j) from (|3.39| ) via a 



generalization of ( p.40[ ) with higher powers of p^ in the arguments of the logarithms. 

However, one quickly realizes that the appearance of a confining phase, i.e. a 
solution of the stationary point equations for 0j of the form of eqs. ( p.48D and ( |3.49| ) 
depends only on the infrared behaviour of the propagators P^^{p'^) for p'^ —* 0. This, 
in turn, is unchanged by the inclusion of operators Oi in ( [5. 61 ) involving higher 
derivatives. 

The general features of the model will thus remain untouched by the inclusion 
of more general operators of the form (|5.5|) , but its "fine-structure" will be affected: 
Lorentz scalar auxiliary fields (pi associated to Lorentz scalar operators Oi involving 
higher derivatives will generally have non-vanishing vevs and affect the vacuum 
structure, and auxiliary fields (pi with Lorentz indices will affect the "response" of 
the model to sources JF,fiu, JB,fiu- 

Since the inclusion of a large number of additional auxiliary fields (pi is evidently 
a rather ambitious program^ the question arises whether - and to what extent - 
one could make use of the property of "universality" of local quantum field theories 
as the ones presented here which, by definition, corresponds to features which are 
practically independent of the UV cutoff A and the bare action. The phenomenon 
of universality arises in the response of a quantum field theory to external sources 
which induce only "low energy" phenomena: The momenta g, of the sources have 
to satisfy qf <^ A^, and the change AE of the vacuum energy (in the case of, e.g., 
constant sources) has to be much smaller than A. (Otherwise one is sensitive to the 
non- universal regime of the effective potential.) 

First, the Wilson loop corresponds to a source which satisfies none of these 

criteria: Even after rewriting it as a source located on a surface S bounded by 

the loop (as in (|4.16|) ) this source varies rapidly in x± perpendicular to S, and its 

^The introduction of pseudoscalar operators of the form e^upcrF^iyPp^ etc., without higher 
derivatives, should however be feasable. We would not expect vevs of the associated auxiliary 
fields, which would spontaneously break parity, but their two point functions would give us infor- 
mations on the topological susceptibility. 
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expectation value depends thus on the non-universal terms involving high derivatives 
in the effective action. Fortunately this does not affect the area law but, as we 
have seen in section 4, the computation of the string tension. In addition, even 
after replacing an infinitely thin surface 5 by a layer of finite thickness, the field 
configurations inside this layer will generically depend on the non-universal regime 
of the effective action unless the "source" is unconventionally weak (or "diluted"). 

The response of the present class of models to sources which vary sufficiently 
slowly in x and which are sufficiently weak will, on the other hand, not depend on 
details of Gbare and/or higher derivative terms. Light quarks in QCD (with extended 
wave functions) could play the role of such sources, and their coupling to the present 
model will constitute an interesting task in the future^. 

To conclude, we have investigated a class of four- dimensional gauge theories with 
finite UV cutoff, which exhibit confinement and allow nevertheless - in the large A^ 
limit - for controllable computations notably in the infrared regime. Some features 
of the confining phase correspond quite to our expectations, notably the possiblity 
to perform a duality transformation of the low energy part of the effective action 
and thus to interpret confinement as monopole condensation. A technically related 
but nevertheless unexpected feature is the appearance of a low energy vector gauge 
symmetry, which allows to "gauge away" the low momentum modes of the abelian 
gluons A^ and explains the absence of the corresponding asymptotic states. If we 
assume in addition, as in eq. ( |2.15| ), higher derivative terms bilinear in 5" such 
that the corresponding massive poles in the propagators are absent, the model has 
no asymptotic states at all (even no bound states, cf. the end of section 5). 

The only "meaning" of the model is thus its reaction to external sources or cur- 
rents. In the case of conserved currents (with respect to the approximate vector 
gauge symmetry) confinement is then actually not based on a long-range attractive 
force between two distinct sources - the corresponding components of the propaga- 
tors cancel precisely for conserved external currents as explained in the paragraphs 
below eq. ( [4.20|) - but on a different phenomenon: The minimal geometric object 
which can be associated to a conserved current or source J^^ is a surface S (cf. the 
expressions (|4.16|) and ([4.19|) ) and, as discussed in section 4, the contribution to the 



action of a source which is constant on 5* is proportional to S (at least for large 

S) inducing the area law of the Wilson loop. The tensor structure T2^^u^pa in the 

^At least a computation of Tf.ff[(j)i, A, B) in a systematic expansion in {dfj_/M-\^), beyond the 
large TV limit, should be feasable. 
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propagator P^^ (cf. eqs. ( |A.2|) , ([A.?]) ) is essential for this result. The question why 
a pointlike particle like a heavy quark can be associated to source for B^y which is 
constant on a surface S is returned to the underlying microscopic theory, such as the 
non-abelian sector of SU{N^ Yang-Mills, but here we found that this is obligatory 
in order to be able to satisfy the constraints ( [4.22| ) 



Interestingly, if the lattice results ||T0[ on the masses of the non-abelian gluons 



in the MAG are correct, an elaboration of the precise relation between the present 
"effective low energy theory" and SU{Nc) Yang Mills - i.e. the structure of S^Jf^ as 
well as the structure of the sources - would involve only physics at comparatively 
small scales of 0{{1 GeV)~^) and is thus perhaps feasable. 

Apart from this evident application of our class of models we would like to point 
out that generalizations in the sense of [|, ^ suggest itself: It is fairly straightfor- 



ward to vary (notably increase) the rank of the tensors A^, B^^, (maintaining rank 
(B) = rank (A) + 1) and to vary the dimension of the space-time. This way one 
can develop dynamical models for the condensation of topological defects of various 
dimensions (p - Branes) in various space-time dimensions. 
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A Appendix 

In this appendix we give explicit expressions for the different propagators, for 
constant fields 0j and in the Landau gauge, both in momentum space and space- 



time and for various limits. For the cutoffs terms ASt (cf. eq. ( |2.1CI| )) we use the 
simplest form 

(A.l) 
(This form of the cutoff on the fields A" is consistent only in the Landau gauge.) 
Subsequently it will be convenient to define the following three 4-index tensors: 



Tl,^u,pa{p) = ^fjipPuPa - Sf^aPuPp - ^z^p PpPa + ^ua PfiPp , 
-'■2,^,1', pa O^pOycr O^crOyp , 

T3^fj,u,pcr{P) = T2^^u,pu 2 Fl,pu,pa\P) ■ (^-2) 

First we give the three propagators P^^, pAB ^^^ pBB ^^ momentum space: 



pA(p2)^^P'+4(03 + i?^(p')) 



Ap^diip^) 
di{p') = crp' (01 + R^ip')) + 4 (01 + R^ip')) U, + R^{p' 



'21 



pZ^p) = ^{Pp5p.-pjpp)p^^{p') 



p^'^ip') 



2p^di{p^) 



p^^ip) = Ti,pu,p. ■ p^^(p') + n,,u,p. 

^^^ Ap^d,{p^) ' 
rf2b') = V + 4i?^(/)+403. 


1 


d2{p') ' 



(A.3) 

For later use it is convenient to introduce the propagators P^^, P^^ defined as 
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PZM = P[,Pu]i p.] = Ti,,uMp) P\p') , (A.4a) 

PZM = PI^PZ^ = Ti,,.Mp) P^^'ip') ■ (A.4b) 

Many physical phenomena induced by the effective action of our model are best 
understood in terms of the propagators in ordinary (Euclidean) space-time in the 
limit where the infrared cutoff is removed. For simplicity we give these propagators 
also for vanishing UV cutoff (which is anyhow non- universal); their short- distance 
singularities should be regularized correspondingly. 

The large distance behaviour of the space-time propagators depends often cru- 
cially on the combination S = (40i03 — (j)l)/a(j)i (defined also in eq. ( p.43D ), which 
vanishes in the "confining phase" (cf. section 3). For completeness we give the 
space-time propagators both for finite S and for S — i> 0. 

First, for finite S, we have 






-^^^^(e€^+^^^i^i-eI^^^^i^i^)) 



+T^^—Tl,^,u,pa{^z) < --^ - —\—Ki \z\2] 






Ti^^u,paidz) is defined in ( |A.2|) with the replacement p^ -^ d^ and Ki is a Bessel 



function. The propagators P^^ and P^^ are 
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^fiu,pa\^) 



^vr^a^i 



-^ l,i.iu,pa\(^z 



Ski 



-i-K,i\z\V^)] . (A.6b) 



Finally, in the limit S ^ the space-time propagators become (we also replace 
^2 by 2 a/0103) 



p;^,^i^) 



1 203 



1 
~2 



^M- I 1^12 



cr 



-d^d^ log kl 



03^1 

2a 



log \z\ + const. 



-(log \z\ + const.') 



P^» 



57r^(TV01 
1 



(c^p^M'^ -dJ^p)\og\z\ 



SttV 



Tl,;.i.,pa(9^)logk| 



1 



167r2( 



"-' l,^u,pa\S^z 



^ ^ ""^-i^J|d2/'^^ 



\z\ V /i 




27r2/i3/2 



Tc 



2,fj.u,pa'i 



-Ki \z\2\ 



(A.7) 



The expressions const, and const.' in ( |A.7| ) actually diverge like log(S). A finite 
result is obtained, however, for P^^ and P^^ in the limit S ^ 0: 



^pV,pCT\^) 



lQ'K^(f)l 



^ T (Fi^) ^ 203, I ,1 

2— T^i,p.,pa(<9,) ^ - ^ log kl !> , 



\z\ 



a 



(A.8a) 



P^^ iz) 



in^'a V 01 



Ti^pu,pa{dz)log\z\ 



(A.8b) 
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